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Abstract 

We present the first global analysis of the Constrained NMSSM that investigates the im- 
pact of the recent discovery of a 126 GeV Higgs-like boson, of the observation of a signal for 
BR (B s — > and of constraints on supersymmetry from ~ 5/fb of data accumulated at the 

LHC, as well as of other relevant constraints from colliders, flavor physics and dark matter. We 
consider three possible cases, assuming in turn that the discovered Higgs boson is: i) the lightest 
Higgs boson of the model; ii) the next-to-lightest Higgs boson; and hi) a combination of both 
roughly degenerate in mass. The likelihood function for the Higgs signal uses signal rates in the 
77 and ZZ — > 41 channels, while that for the Higgs exclusion limits assumes decay through the 
77, tt, ZZ and W + W~ channels. In all cases considered we identify the 68% and 95% credible 
posterior probability regions in a Bayesian approach. We find that, when the constraints are 
applied with their respective uncertainties, the first case shows strong CMSSM-like behavior, 
with the stau coannihilation region featuring highest posterior probability, the best-fit point, a 
correct mass of the lightest Higgs boson and the lighter stop mass in the ballpark of 1 TeV. We 
also expose in this region a linear relationship between the trilinear couplings of the stau and 
the stop, with both of them being strongly negative as enforced by the Higgs mass and the relic 
density, which outside of the stau coannihilation region show some tension. The second and the 
third case, on the other hand, while allowed are disfavored by the constraints from direct de- 
tection of dark matter and from BR (B s — > Without the anomalous magnetic moment 
of the muon the fit improves considerably, especially for negative effective fi parameter. We 
discuss how the considered scenarios could be tested furtherly at the LHC and in dark matter 
searches. 



*On leave of absence from the University of Sheffield, UK. 
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1 Introduction 



In July 2012 the CMS and ATLAS collaborations at the LHC made the announcement of a ~ 5a 
discovery of a boson with mass 125.3 ± 0.6 GeV [1] and 126.0 ± 0.6 GeV [2], respectively, consistent 
with the Higgs boson predicted by the electroweak (EW) Standard Model (SM). The CMS Col- 
laboration has recently updated its results [3] to 125.8 ± 0.6 GeV. The updated result is based on 
the analysis of the data corresponding to integrated luminosities of 5.1/fb at ^/s = 7TeV and up 
to 12.2/fb at yfs = 8TeV in the 77, ZZ, WW, tt and bb decay channels. The ATLAS analysis 
combined approximately 4.8/fb of 2011 data at yfs = 7TeV in the same five channels with 5.8/fb of 
data at y/s = 8 TeV in the ZZ, 77 and WW channels only. Evidently, the excess of events near the 
reported mass is driven by the 77 and ZZ channels, owing to the highest mass resolution of these 
channels. Both collaborations also determined the ratio fJ,(X) of the observed Higgs production 
cross section to the one predicted by the SM, in each of the mentioned Higgs decay channels X. 
An enhancement in the 77 channel was reported by ATLAS, with ^(77) = 1.9 =L 0.5, as well as by 
CMS, with ^(77) = 1.6 ± 0.4. On the other hand, the updated value of (J,(ZZ) by CMS @] and 
ATLAS [5] is, within la error, SM-like. 

While its exact characteristics are still being carefully examined, such a boson is not only con- 
sistent with the SM Higgs particle but can also be easily accommodated into models of new physics, 
particularly those based on softly-broken low-energy supersymmetry (SUSY) which actually pre- 
dict a relatively light Higgs boson. Most studies have been performed within the framework of 
the two Higgs doublet Minimal Supersymmetric Standard Model (MSSM), usually augmented by 
various unification boundary conditions. The framework, however, suffers from the so-called '/u- 
problem' [6] : the SUSY-preserving parameter [i in the superpotential is expected to be of the order 
of the SUSY breaking scale Msusy ~ 1 TeV to assure correct radiative EW symmetry breaking. 

Perhaps the most compelling and simplest solution to the /i-problem is to invoke an additional 
gauge singlet field coupled to the Higgs doublets of the MSSM [6|. The //-term is then generated 
dynamically through the vacuum expectation value (vev) of the singlet field which is of the order 
of Msusy- The model is commonly referred to as the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM); for reviews see, e.g., 0[B]. After the discovery of the Higgs-like boson, numerous 
studies have appeared in the context of the NMSSM [9j [TOj, [Til G21 [13], since the model presents 
several interesting features. In the NMSSM the particle content remains the same as in the MSSM, 
except that the number of Higgs bosons increases from three to five: three neutral scalars ^1,2,3 and 
two neutral pseudoscalars a\^- The number of neutralinos also increases from four to five, Xi,...,5> 
owing to the singlino partner of the singlet Higgs boson. 

On the phenomenological front, a two-loop theoretical upper bound on the mass of the lightest 
CP-even Higgs boson, m/ 11 , can increase by a few GeV compared to the MSSM for some com- 
binations of the model's parameters. Moreover, in the context of the recent LHC results, the 
next-to-lightest CP-even Higgs boson of the model, /12, could also have mass around 126 GeV 
with h\ being even lighter than the LEP bound on SM-like Higgs boson mass, without violating 
it [THUS], owing to the singlet-doublet mixing effects. Both features potentially allow to achieve 
the correct mass of the experimentally detected Higgs without excessive fine-tuning [To ] fT7] fT8 } [TP] . 
making the NMSSM potentially a more 'natural' model than the MSSM. Finally, the NMSSM also 
offers the additional possibility that the observed excess in the 77 and ZZ rates could be due to the 
combination of h\ and /12 decays (with and mh 2 being almost degenerate) [ID]. Besides masses, 
the decay rates of the Higgs boson(s) are also affected by the modifications in the superpotential. 
For example, for large doublet-singlet mixing, the h\bb coupling becomes suppressed, reducing the 
total decay width. As a result, the branching ratio (BR) of h\ — > 77 becomes marginally enhanced 
compared to the SM or the MSSM [20]. An important consequence of this feature is that hi lighter 
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than 115 GeV could have escaped detection at LEP, and also at the LHC, due to the possibility 
of its decay into a pair of the two lighter pseudoscalars, hi — > a\a\ (see, e.g., [21]), causing thus a 
suppression in all of the other BR's of h\. 

The NMSSM in its most general form contains more than a hundred free parameters: those 
appearing in the MSSM and some additional ones relating to the extended Higgs sector. Similarly 
to the MSSM, where imposing universality conditions on the soft SUSY-breaking parameters of the 
model at the scale of grand unification (GUT scale) leads to the Constrained MSSM (CMSSM) [22J, 
a GUT-constrained version of the NMSSM (CNMSSM) [231 [Ml US] can be defined; see Sec. for 
more details and discussion. 

It is well known that, in the CMSSM the lightest Higgs boson's mass, as calculated at 2-loop 
level with FeynHiggs [26] or SoftSusy [27], is typically a few GeV below 126 GeV, which can be 
considered as somewhat unsatisfactory even if one takes into account a residual error of some 2- 
3 GeV [28] from scheme dependence. (Values of nih around 126 GeV can still be obtained in the 
CMSSM for Msusy ^ 1 TeV at the expense of increased fine-tuning [29].) A priori one could 
expect that in the CNMSSM, with more freedom in the Higgs sector, the tension will be reduced. 
However, recent studies [30] using a fixed-grid scan approach have shown that in the CNMSSM it is 
extremely difficult to obtain a Higgs boson as heavy as 126 GeV, particularly one with an enhanced 
77 decay rate, while also satisfying other phenomenological constraints, thus nullifying the noted 
advantages of the singlet-extension of the MSSM. 

On the other hand, it has been demonstrated by our recent global analysis of the CMSSM |29j 
and by some earlier Bayesian studies \31\ 132] , that a proper treatment of the experimental con- 
straints through a likelihood function can lead to significantly different results from fixed-grid 
scans where such constraints are typically implemented in a more simplified box-like fashion, with 
all points accepted when satisfying experimental values within some fixed range (typically within 
la, or so), and otherwise rejected. One of the main advantages of the statistical approach is that 
scanned points are instead 'weighted' by the total x 2 > thus indicating how well they fit all con- 
straints. For example, in the fixed-grid approach a point giving a value of even one constraint 
slightly beyond the allowed range while reproducing central values of all the other constraints 
would be rejected, while a point with values for all the constraints barely within the allowed boxes 
would be completely allowed. In contrast, in a statistical approach both points would be accepted 
but weighted with their respected \ 2 - Another advantage of the statistical approach based on the 
likelihood function is that theoretical and experimental uncertainties can be easily implemented in 
a consistent manner. Thanks to some recent developments in sampling algorithms (see, e.g., |33] 
and |34|). multidimensional Bayesian scans can now be carried out rather quickly and efficiently. 

One should note two important changes in the data that have recently taken place on the 
experimental front. In November 2012 the LHCb Collaboration published the most recent update 
of their search for the rare decay B s -4 [35], reporting an excess of B s — > candidates 

over the background. The measured value, BR(B S — > = 3.2^2 x 10~ 9 , is now very close 

to the time-averaged SM value, 3.5 x 10 -9 |36l I37| . We include the constraint in our likelihood 
function taking into account both theoretical and experimental uncertainties, as will be described 
below. 

The other important update was the top pole mass by the Particle Data Group, obtained from 
an average of data from Tevatron and the LHC at yfs = 7 TeV, Mt = 173.5 ± 1.0 GeV |38j . As we 
shall see below this is a welcome increase relative to its previous value in the context of the Higgs 
sector of constrained SUSY models as it pushes the mass of h\ up, closer to the experimentally 
observed Higgs- like resonance mass. 

In this article, we present the first global Bayesian analysis of the CNMSSM after the observation 
of the SM Higgs- like boson. We separately consider the cases of this boson being h\, or /12, or a 
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combination of both. We test the parameter space of the model against the currently published, 
already stringent constraints from SUSY searches at the LHC and other relevant constraints from 
colliders, 6-physics and dark matter (DM) relic density. Our goal is to map out the regions of 
the parameter space of the CNMSSM that are favoured by these constraints. As in our CMSSM 
study [29], the CMS razor limit based on 4.4/fb of data is implemented through an approximate 
but accurate likelihood function. We also study the effects of relaxing the (g — 2)^ constraint. 

The article is organized as follows. In Sec. [2] we briefly revisit the model, highlighting some of 
its salient features. In Sec. [3] we detail our methodology, including our statistical approach and our 
construction of the likelihoods for the BR (B s — > signal, the CMS razor 4.4/fb analysis, and 

the CMS Higgs searches. In Sec. [4] we present the results from our scans and discuss their novel 
features. We summarize our findings in Sec. [5} 

2 The NMSSM with GUT-scale universality 

The NMSSM is an economical extension of the MSSM, in which one adds a gauge-singlet superfield 
S whose scalar component couples only to the two MSSM Higgs doublets H u and H d at the tree- 
levelQ The scale-invariant superpotential of the model has the form 

W = XSH u H d + *-S 3 + (MSSM Yukawa terms) , (1) 

where A and k are dimensionless couplings. Upon spontaneous symmetry breaking, the scalar Higgs 
field S develops a vev, s = (S), and the first term in Eq. assumes the role of the effective /i-term 
of the MSSM, ^ e g = As. The soft SUSY-breaking terms in the Higgs sector are then given by 

V so{t = m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 \S\ 2 + {\A x SH u H d + ^ K S 3 + h.c.J , (2) 

where A\ and A K are soft trilinear terms associated with the A and k terms in the superpotential. 
The vev s, determined by the minimization conditions of the Higgs potential, is effectively induced 
by the SUSY-breaking terms in Eq. ([2]), and is naturally set by MgusYj thus solving the //-problem 
of the MSSM. 

The Constrained version of the NMSSM (CNMSSM) is defined in terms of five continuous input 
parameters and one sign, 

mo,m 1/2 , A),tan/3, A,sgn(^ efr ) , (3) 

where unification conditions at a high scale require that all the scalar soft SUSY-breaking masses 
in the superpotential (except mg) are unified to tuq, the gaugino masses are unified to rra^, and 
all trilinear couplings, including A\ and A K , are unified to Aq. Additionally, tan/3 (the ratio of 
the vevs of the neutral components of the H u and H d fields) and sgn(^ e fj) are taken as input at 
the weak scale, and A is an input parameter at Msusy- The singlet soft-breaking mass m 2 s is not 
unified to uiq for both theoretical and phenomenological reasons. From the theoretical point of 
view, it has been argued [39J that the mechanism for SUSY breaking might treat the singlet field 
differently from the other superfields. From the phenomenological point of view, the freedom in rag 
allows for easier convergence when the renormalization group equations (RGE's) are evolved from 
the GUT scale down to Msusy- It also yields, in the limit A — > 0, and with As fixed, effectively the 
CMSSM plus a singlet and singlino fields that both decouple from the rest of the spectrum. 

For simplicity we will be using the same notation for superfields and their bosonic components. 
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See text 
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0.0056, 10% 
8.0, 1.0 
0.22, 0.21 
0.66, 0.38 

0.043 ps" 1 , 2.400 ps" 1 
+1.5- 1.2, 10% 


Gaussian 
Gaussian 
Gaussian 
Gaussian 
Gaussian 
Gaussian 


ED 

[46] 

S3 
EH] 
[35] 



Table 1: The experimental constraints that we apply to constrain model parameters. mh aig , m h hi 
R heig and R hh 



id ' 



are defined in Sec. 3.3 



3 Statistical treatment of experimental data 

We explore the parameter space of the model with the help of Bayesian formalism. We follow the 
procedure outlined in detail in our previous papers |40[ HI] [29], of which we give a short summary 
here. Our aim is to map out the 68% and 95% credible regions of p(m\d), the posterior probability 
density function (pdf), given by Bayes' theorem, 

p(d\£(m)) = C is the likelihood function, which describes the probability of obtaining the data d 
given the computed value of some observable £(m), which is a function of the model's parameters m. 
C also incorporates the experimental and theoretical uncertainties. Prior probability 7r(m) encodes 
assumed range and distribution of m. Finally, p(d) is the evidence, which is a normalization 
constant as long as only one model is considered, but serves as a comparative measure for different 
models or scenarios. 

Bayes' theorem provides an efficient and natural procedure for drawing inferences on a subset 
of r specific model parameters (including nuisance parameters), or observables, or a combination 
of both, which we collectively denote by ipi. They can be obtained through marginalization of the 
full posterior pdf, carried out as 

P^i=i,..,r\d) = J p(m\d)d n - r m, (5) 

where n is the total number of input parameters. An analogous procedure can be performed with 
the observables and with a combination of the model's parameters and observables. 

In order to evaluate the posterior probability given by Eq. Q, one needs to first construct the 
likelihood function. The constraints that we include in the current analysis are listed in Table [TJ 
We shall be discussing them in turn below. As a rule, following the procedure developed earlier [48J, 
we implement positive measurements through the usual Gaussian likelihood, while upper or lower 
limits through an error function smeared with both theory and, when available, experimental error. 
The construction of the likelihoods for direct SUSY and Higgs searches is more involved, and will 
be explained in detail later in this section. 



5 



3.1 Likelihood for BR (B s -> fj+fj, ~) 

In November 2012 the LHCb Collaboration published the most recent update of their search for the 
rare decay B s — > fi + /i~ [S3], based on a combination of the 2012 data samples of 1.1/fb of proton- 
proton collisions at \/s = 8TeV and the 2011, 1.0/fb data at \fs = 7TeV. The data superseded 
the combination of 2011 data from ATLAS, CMS and LHCb published in June 2012 [35] . 

LHCb observed an excess of B s — > candidates over the background, consistent with the 

SM expectation. The measured value is BR(B S — > = 3.2jl^'| x 10 -9 , with a statistical 

significance of 3.5 a. We used this information to construct an approximate likelihood function 
for BR (B s —7- which parametrizes the experimental and theoretical uncertainties. The most 

important sources of theoretical uncertainty are the B s decay constant (main contribution) and 
its lifetime, the top quark mass and the CKM matrix elements Vj* b Vt s [37J and their total amounts 
to approximately 11% of the mean value [50]. Thus, we parametrized the likelihood function 
as a combination of two half-Gaussians, to take into account the asymmetry in the experimental 
uncertainty. A theoretical uncertainty equal to 10% of the calculated value was added in quadrature. 
Notice that we neglected the uncertainty due to the top pole mass (~ 1%) since it is taken care of 
by scanning over the SM nuisance parameters (see below). 

3.2 Limits on SUSY from the LHC 

In Ref. [29] we derived a methodology for constructing an approximate likelihood map that repro- 
duced the lower limit in the (mo, m 1 ^ 2 ) plane set by the CMS Collaboration by applying the razor 
method [511 [42] to SUSY searches in all-hadronic modes based on 4.4/fb of data. We validated 
our map against the experimental results by evaluating the resulting 95% CL contour in the (mo, 
m 1( / 2 ) plane and comparing it with the corresponding 95% CL exclusion limit provided by the CMS 
Collaboration. We obtained a very good agreement confirming that our procedure for generating 
these likelihood maps was indeed correct. Our methodology can be applied to produce the SUSY 
exclusion limits in any .R-parity conserving supersymmetric model, as long as the supersymmet- 
ric spectra in that model present similar features to the ones of the CMSSM, namely a sufficient 
mass difference between the lightest SUSY particle (LSP) and gluinos or squarks of the first two 
generations. When extending the procedure to other models, one should also take into account 
possible changes in the production and decay modes of the particles. In this paper we follow the 
same methodology for generating a SUSY likelihood map based on the CMS razor 4.4/fb analysis 
for the NMSSM. 

For each point in an mo — mi/2 grid, with 50GeV step size in both parameters, the super- 
symmetric mass spectrum and decay table are generated using NMSSMTools [52], and fed into 
Herwig+- 1- [53] for parton shower generation and calculation of the cross sections. Herwig++ al- 
lows one to work in the framework of the NMSSM, and hence takes care of a possible contribution 
from the extended neutralino and Higgs sectors. The production modes are not influenced by these 
sectors, since the razor search is designed to detect squarks and gluinos, which are produced through 
color interactions. The output of Herwig-(— |- is passed in HepMC [53] format to Delphes |55j for 
reconstruction of the physical objects by simulating the CMS detector's response. The output from 
Delphes is distributed into 38 bins, defined in terms of the razor variables Mr and R 2 , each con- 
taining the data (a number of background and signal events) from the experiment, to produce an 
efficiency map for each bin, which is then translated into a likelihood map for the whole grid. We 
refer the reader to Ref. [29] for more details of the complete procedure adopted for the production 
of efficiency maps. 

The 95% CL razor contour thus reproduced is shown in Fig. [T] The actual 95% CL line obtained 
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Figure 1: The 95% CL lower bounds from our approximate razor likelihood, for different values of 
A and A K , compared with the experimental line (in dashed black). 



by CMS is also shown as a dashed black line. The relevant CNMSSM parameters were fixed as: 
tan/? = 5, Aq = 0, A K = —15 GeV, A = 0.001 and sgn(^ e g) = +1. The produced line is indicated 
in purple. One can see that the CMS limit is reproduced very well in the limits of m^rt 3> ttiq 
and vice versa. A slight disrepancy in between those two extremes is basically irrelevant, as we 
shall see below. The dependence of the limit on A and A K is negligible, since these parameters only 
affect the boundaries of the physically allowed region but do not cause a shift in the position of 
the contour. This is clear from Fig. [I] where we show a contour for A = 0.1 (in green) and another 
for A K = —1000 GeV (in red), fixing the other parameters to the stated value in each case. The 
negligible dependence of the razor limit on Aq, tan/3 and sgn(/x e ff) has already been verified in 
many analyses; see e.g., (56] [571 [29] . Our choice of these parameters thus only ensures maximum 
allowed physical region in the (mo, m 1 / 2 ) space, since the physicality aspect is later taken into 
account during the scan when the unphysical points are discarded. 

3.3 Likelihood for the Higgs bosons 

The measured mass of the Higgs-like boson has important consequences for global properties of the 
CNMSSM. 

The upper limit on the mass of hi is given by [58] 

\2 

m 2 hl < Ml cos 2 2(3 + A V sin 2 2/3 - -^v 2 

where Am 2 denotes loop corrections, and v = 174 GeV is the tree-level Higgs doublet vev. The 
second, NMSSM-induced term on the right-hand side of Eq.Q can enhance the hi mass by a few 
GeV beyond its MSSM value, if A is large enough, A > 0.5, and tan/3 < 2. For larger values of 
tan/3, the second term on the right-hand side becomes suppressed, and the third, negative term 
becomes significant. Then, the value of the Higgs boson mass can only be enhanced by the loop 
corrections. As we shall see, when other constraints are taken into account, especially the ones 



A — sin 2/3 ( k 



Ax 
2s 



+ Am 2 , 



(6) 
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coming from the LEP Higgs searches, 6-physics, the muon anomalous magnetic moment, {g — 2)^, 
and DM relic density, the aforementioned conditions on A and tan (3 become hard to satisfy and 
drastically restrict the parameter space of the model. 

As noted in the Introduction, in the NMSSM either of the two lightest CP-even Higgses could 
serve the role of the SM-Higgs-like boson observed at the LHC. There are thus three interesting 
possibilities entailing such a discovery^] 



Case 1 
Case 2 
Case 3 



m hl ~ 125.8 GeV, 
m h2 ~ 125.8 GeV, 
m/jj ~ rrifo ~ 125.8 GeV. 



3.3.1 Mass Likelihood 

In this paper we analyze all the three cases mentioned above. We will henceforth refer to our 
selected 125.8 GeV signal Higgs for either of the first two cases generically as h s i g . We impose a mass 
requirement on h S i g through a Gaussian likelihood which takes into account both the theoretical 
(r) and experimental (a) errors, following the procedure detailed in Ref. [29J, 

£ffig g s ma^m^J = exp [-(125.8 GeV - m hsis ) 2 /2(r 2 + a 2 )] . (7) 

The theoretical error is due to residual differences between calculations using different approaches 
and schemes, and it is estimated in the literature to be of the order of 2 — 3 GeV [59J. We assume 
r = 3 GeV. We use the CMS determination of the Higgs mass rather than the one by ATLAS, 
since it has been most recently updated, and also for consistency with our previous work [29J. 



3.3.2 Cross section Likelihood 

In order to be thoroughly consistent with the CMS measurement, we need to calculate for both 
Higgs bosons %,2 the reduced cross sections, defined in the literature as 

p /yx <j{pp^hj) BR(hj^X) 

a{pp -> /ism) BR(h S M -> X) 

for a given Higgs decay channel, X. 

Equation ^ can be approximated by 

R hi (X)= Rl(x)n SM (Y), (9) 

Yd prod 

where the sum runs over the Higgs production channels Y (Y = gg for gluon-fusion, VV for vector 
boson-fusion and Higgs-strahlung off a Z boson, tt and bb for associated Higgs production with top 
and bottom quarks, respectively), and the ratios Hsm(Y) = o{j>p — )■ Y — > ^sm)/c(pp — > ^sm) are 
obtained from the public tables provided by the LHC Higgs Cross Section Working Group [SH |UT] 
for = 8TeV. 

The reduced cross sections RX (X) for the individual production channels are provided by the 
mass spectrum generator, NMSSMTools, that we use for our analysis. Alternatively, they can be 
expressed in terms of the Higgs reduced couplings C(X) (the ratio of the couplings of the Higgs 
boson with a given mass to a pair of X particles in the NMSSM, to the ones calculated in the 



2 We assume that the third CP-even Higgs is always heavier than the current experimental reach. 
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SM) and the decay branching ratios BR(/i 
generator, 



Rh. 



(X) 



X), which are also provided by the mass spectrum 
a(Y -> hi) BR(hi -> X) 



C 2 (Y) x 



^sm) BR(/ism — > 

r(/ li x)/r tot 



r(/t 



SM 



x)/r SM 



c 2 (y)c 2 (x) 



Fe SM decay 



tot 

BR(hi 



F) 



C 2 (F) 



where the sum runs over the decay channels F open to the SM Higgs boson. 

We require Rh Bi (77) and Rh si (ZZ) to comply with the measured rates, i? Q bs(77) 



(10) 



1.6±0.4 [I] 



and R ohs (ZZ) = 0.8±g;|| @J, respectively^ For the 77 channel the 'signal' likelihood is taken to be 



a Gaussian around the measured central value, 



£si g (77) = ex P ~ (^obs(77) - Rh Big {ll)) /2(ct 2 7 + r 2 7 ) 



(11) 



with o" 77 being the experimental error. For the ZZ channel, since the experimental error is asym- 
metric, the signal likelihood is defined in terms of half-Gaussians for the positive error a% z and neg- 
ative error o~^ z each, as done in the case of the B s — > likelihood above, tx ~ 15% x R^. (X), 
is a very conservative estimate of the theoretical error based on [62], used for every channel X. 

In addition to constraining /i s ; g , there is another crucial aspect of the Cases 1 and 2, which is 
that the second of the two lightest CP-even Higgs bosons must have escaped detection at the LHC, 
or be 'hidden'. In the following we refer to it as /ihid- In other words, the production rate of /ihid 
should be less than what the experiments are currently sensitive to for all X. As a result, Rh hid (X) 
should also be constrained by experimental data, so that all points where the rate of /ihid is large 
enough for it to have been observed, are rejected in our analysis. For this purpose, we construct an 
'exclusion' likelihood. Following the procedure outlined in |48| for exclusion bounds we first define 
a step function, 



^excl ( m K^ R h h 



AX), m (X)) 



(12) 



1 for R hbU (X) < m {X) 
for R hhlA (X) >// 95 P0, 

where ^95 (X) is the value of the signal strength modifier fJ>(X) = (Jh bid (X) /<Tft SM (X) that is excluded 
at 95% CL by the LHC searches for a Higgs with a given mass run = m hud> obtained from the 
latest exclusion plots published by the CMS collaboration [63] j^] Then, in order to take into account 
the theoretical error on the true values of the reduced cross section and the Higgs mass, £ 



(step) 
excl 



IS 



smeared out further by convolving it with Gaussian functions centered around their true theoretical 
values Rh hid (X) and rhh hid , respectively, so that the exclusion likelihood now becomes 



C 



(smear) 
excl 



RhwV) = I dm hh . d I dR hidd C^\m hhid ,R hhid ,fi) 



x exp 



(m hhid - m hhic 
2r 2 



\2l 



exp 



(Rh 



Rh h 



2f 2 



(13) 



We only use the dominant decay channels where a ~ 5a excess has been observed. 
4 We used the latest exclusion limits provided based on 5.1/fb of data at yfs — 7TeV and 12.2/fb of data at 
y/I = 8TeV for the ZZ channel, on 4.9/fb at = 7TeV and 12.1/fb at -/a = 8TeV for the WW -> llw channel 
and on 17/fb at y/s — 7 — 8TeV for the rr channel. Limits for the 77 channel have not been updated by CMS and 
are, therefore, still based on a combination of 5.1/fb at y/s = 7TeV and 5.3/fb at y/s = 8TeV of data. 
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± 11U1 lylBlllULLHUll 




TTnivprQftl Qpalar Triple 

111 V L I OCLl Ot-jCLlCtl lllcLoo 


inn 4000 




m l/2 


Universal gaugino mass 


100, 2000 


Log 


A 


Universal trilinear coupling 


-7000, 7000 


Linear 


tan /3 


Ratio of Higgs vevs 


1, 62 


Linear 


A 


Higgs trilinear coupling 


0.001, 0.7 


Linear 


Nuisance 


Description 


Central value ± std. dev. 


Prior Distribution 


M t 


Top quark pole mass 


173.5 ± 1.0 


Gaussian 


m b (m b ) MS 


Bottom quark mass 


4.18 ±0.03 


Gaussian 


as(M z ) MS 


Strong coupling 


0.1184 ± 0.0007 


Gaussian 



Table 2: Priors for the parameters of the model and for the SM nuisance parameters used in our 
scans. Masses and Aq are in GeV. 



where the theoretical errors are taken to be r = 3 GeV and f = 15% • Rh hid [62] , respectively. The 
exclusion likelihood is calculated for X = 77, ZZ, WW and rr. Finally, in order for our exclusion 
criterion to be consistent with our criterion for signal observation at 125.8 ± 3.1 GeV (with theory 
and experimental errors added in quadrature), we further impose the condition 

J for 122.7 GeV < m hhid < 128.9 GeV 

£exdW - j elsewhere. (M) 

For Case 3, the likelihood in Eq. Q is separately computed for both hi and /i2- In this 
scenario, only the combined production rate for hi and /12 needs to be equal to R \>s(X). Hence 
the observation likelihood is now defined as 

C ohs (X) = exp{- [R obs (X) - (R hl (X) + R h2 (X))f /2{a\ ±4)} . (15) 
The values and uncertainties of our Higgs constraints are given in Table [TJ 



4 Methodology and Results 

The scanned ranges of the CNMSSM parameters along with the type of distribution of their prior 
are listed in Table [2] Also listed there are the input ranges of the nuisance parameters included in 
the scans. The sign of /i e s is fixed for a given scan and equal to ±1 or —1. 

The analysis was performed using the package BayesFITS which engages several external, pub- 
licly available tools: for sampling it uses MultiNest |64j with 4000 living points, evidence tolerance 
factor set to 0.5, and sampling efficiency equal to 0.8. Mass spectrum, along with AMb s , is com- 
puted with NMSSMTools v3.2.1 and passed via SUSY Les Houches Accord format to Superlso 
v3.3 [65| to calculate BR (B X s7 ), BR (B s fi+fj,-), BR (B u tu), and 8 (g - 2)J USY . DM 
observables, such as the relic density and direct detection cross sections, are calculated with Mi- 
crOMEGAs 2.4.5 [66] • 

Below we will present the results of our scans as one-dimensional (ID) or two-dimensional (2D) 
marginalized posterior pdf maps of parameters and observables. For example, in evaluating a 
posterior pdf for a given parameter, we marginalize over the model's all other parameters and the 
SM nuisance parameters, as described in detail in Refs. [29], |4"1]. 

Notice that when discussing the results of the global scan for Case 1 it will become apparent 
that this case presents a remarkable CMSSM-like behavior. It would therefore be natural to try to 
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compare those results with our recent CMSSM analysis [29J. In doing so, one needs to take into 
account the differences between the numerical codes and constraints adopted in both studies. We 
summarize them here. 

1. In the present study we use NMSSMTools for calculating the supersymmetric spectrum, while 
in [29] we used SoftSUSY. We have repeatedly cross-checked the spectra obtained in the MSSM 
limit of the NMSSM with the ones generated by SoftSUSY, finding some differences, especially with 
respect to loop corrections giving the largest values of the lightest Higgs mass. In some regions of 
the parameter space the difference between the two generators can amount to a maximum of ~ 0.5- 
1 GeVj^] Given the experimental and theoretical uncertainties in the Higgs mass, such difference 
translates into ~ 0.25 units of x 2 > which is not significant for the purpose of the global scan. 

2. In this paper we use the value of BR(B S — >• measured at LHCb [35], which has 
been incorporated in the likelihood as described is Sec. |3.1[ The SM rate rescaled by the time 
dependent asymmetries is now BR (B s — > /*i + /OsM = (3.53 ± 0.38) x 10~ 9 [50], which is a value 
more appropriate for comparison with the experimental rate than the unsealed, ~ 3.2 x 10 -9 , one. 

3. We have updated the nuisance parameters Mt and mb(mb) MS following [38]; see Table ij 
The upgrade in Mt has significant implications for m^- The leading one- loop corrections to the 
Higgs mass squared are given by 



Ami - 



4tt 2 v 2 



2 USY \ + _ 

™? J m svsy V 12Mf USY 



i n / Mhm.) + ( i _ 



(16) 



where m; is the running top quark massj^] Mstjsy is the geometrical average of the physical stop 
masses, Msusy = ^ m t 1 m t 2 > an< ^ X t = A t — /j, c g cot /3. Since Am^ oc m\ it is now easier to generate 
Higgs masses in agreement with the experimental values. In particular, as we highlighted in [29], 
a Higgs mass compatible with the observed excess at 126 GeV was rather difficult to achieve over 
the CMSSM parameter space. That tension has now become somewhat reduced, and we will show 
below that the correct Higgs mass can be obtained in the CMSSM limit of the CNMSSM. 



4.1 Impact of the relic density 

To set the ground for the presentation of our numerical results, we first comment on the role of 
the relic density of DM in selecting favored regions. The relic density is a strong constraint, since 
it is a positive measurement (in contrast to a limit) with a rather small experimental uncertainty 
(Table [TJ. On top of it, it is well known that in unified SUSY models with neutralino LSP the 
corresponding abundance Q x h * s typically too large, or in other words, its annihilation in the 
early Universe is 'generically' too inefficient. Specific mechanisms for enhancing it are therefore 
needed which, however, are only applicable in specific SUSY configurations. As a result, in most 
cases the regions of high probability in the global posterior will reflect one or more of the regions 
of parameter space where Q x h 2 is close to the measured relic density of DM. The regions that are 
still allowed by direct SUSY searches are: 

1. The stau-coannihilation (SC) region [68]. As is known, in constrained SUSY models, like 
the C(N)MSSM, this is a narrow strip at a sharp angle to the axis. The values of Aq and 

tan/3 are also constrained, as only for \Aq\ not exceeding ~2-3TeVand tan/3 not too large does 
the mass of the stau become light enough to be comparable with the neutralino mass, but not 
too light so to make it the LSP. Values of m^ 2 that are excessively large, on the other hand, 

5 The best agreement between SoftSUSY and NMSSMTools is obtained by setting the flag precision for the 
Higgs masses to zero which, therefore, was chosen as default setting for our calculations. 

6 Note that running top quark mass is related to the pole mass through the formula given in Eq. (10) of Ref. [67] . 
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Figure 2: Marginalized 2D posterior pdf in (a) the (mo, rn 1 / 2 ) plane and (b) the (Aq, tan/3) plane 
of the CNMSSM for Case 1, constrained by the experiments listed in Table [TJ The 68% credible 
regions are shown in dark blue, and the 95% credible regions in light blue. The solid black (dashed 
gray) line shows the CMS razor (ATLAS hadronic) 95% CL exclusion bound. 



can suppress the annihilation cross section [69] . After other relevant constraints are included, the 
parameters of interest are, therefore, mo < 600 GeV, m 1( / 2 < 1000 GeV, \Aq\ < 3000 GeV and, 
when the neutralino is close to 100% bino, tan/3 < 30. A similar effect can also be obtained for 
large Aq with the stop t\ replacing f\ [70] . 

2. The A-funnel (AF) region, where neutralinos annihilate through the resonance with the 
lightest pseudoscalar |71j . This mechanism can occur over broad ranges of the (mo, mi/ 2 ) plane 
where the pseudoscalar mass is close to twice the neutralino LSP mass, and is enhanced by large 
tan/3 (tan/3 > 35) and positive Aq. 

3. The focus point/hyperbolic branch (FP/HB) region [721 173] . where the annihilation cross 
section can be enhanced by an increased higgsino component of the neutralino. For this to occur, [i 
(or fi e ff in the NMSSM) must be of the order of a few hundred GeV, and tan /3 cannot be too large, 
tan/3 < 45. In the (mo, m 1 / 2 ) plane the condition corresponds to the region where mo 3> rn l / 2 . 

4.2 Impact of the Higgs mass 

The measurement of the Higgs mass has added an important additional constraint on unified SUSY 
models. Below we will be discuss in turn the three cases listed earlier in Sec. 13.31 



Case 1. In Fig. t a) we show the marginalized posterior pdf in the (mo, m^ 2 ) plane of the 
CNMSSM for Case 1, obtained by imposing simultaneously all the constraints shown in Table [TJ In 
these and the following plots the Bayesian 68% (la) credible regions are indicated in dark blue and 
the 95% (2a) credible regions in light blue. Notice that the regions of high probability are located 
above the CMS razor limit, which we implemented in the likelihoood as described in Sec. 3.2, and 
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Figure 3: Marginalized 2D posterior pdf in (a) the (m^, m f - ) plane and (b) the (m^, X t ) plane, 
for the CNMSSM constrained by the experiments listed in Table [T] for Case 1. The color code is 
the same as in Fig. [2j 



which is shown in the plots as a solid black line. In Case 1 the role of the SM-like Higgs is played 
by the lightest CP-even scalar (almost purely .ff^-like), while hi (almost purely i^-like) and /13 
(almost purely singlet-like) are usually much heavier and decoupled. This case is thus expected to 
present features very similar to the CMSSM. 

In Fig. 2;[a) one can see two main 68% credibility regions: the SC region on the lower left 



side, and the AF region on the top part of the plot. As is also the case in the CMSSM, besides 
giving the correct relic abundance, the SC region shows also the better fit to the Higgs mass, 
m/ ll = 124.5 GeV. This is because, as we explained in Ref. [29], in the SC region Aq can easily be 
negative without spoiling the relic abundance constraint. Large negative values of Aq are necessary 
to drive the parameter At to even larger negative values at the EW scale, thus making the stop 
mixing contribution to the loop corrections of the Higgs mass maximal. Note as well that the 
best-fit point is also located in the SC regionj^] In Fig. c ± b) we show the marginalized posterior in 
the (^4o> tan/3) plane. The high probability 'island' at negative Aq and tan/3 < 25 corresponds to 
the SC region. 

In the CNMSSM the SC region appears to be more extended relative to the CMSSM [29], and 
somewhat larger Higgs masses also seem allowed, as we will show below. The increased relevance 
of the SC region is due to the fact that it is now much easier to obtain values of the Higgs mass 
closer to the correct one. This could be mistakenly thought to be a specific feature of the CNMSSM 
extended Higgs sector. We have checked that this is not the case: the Higgs mass is simply quite 
sensitive to the increased central value of the top mass, as we explained at the beginning of this 
section. 

In Figs. 5|[a) and |3|[b)[ we show how the constraints affect the main observables responsible for 



7 We postpone further discussion of the best-fit points until Sec. 
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Figure 4: Marginalized 2D posterior pdf in (a) the (A t , /i e ff) plane and (b) the (At, A T ) plane, for 
the CNMSSM constrained by the experiments listed in Table [T] for Case 1. The color code is the 
same as in Fig. [2j 



the loop corrections to the Higgs boson mass. In Fig. 3;[a) we show the posterior in the (m/ ll 



Trif ) plane. Notice that the lightest stop does not need to be excessively heavy in the SC region, 



m 



< 1 TeV, since large stop mixing compensates for smaller Mstjsy- The correct Higgs mass is 
obtained since \Xt\ ~ 2.5 TeV, as can be seen in Fig. 3;[b)| where we plot the posterior in the (m^, 
Xt) plane. On the other hand, one can see that in both figures the best-fit point lies outside of 
the 68% credibility regions of the marginalized 2D pdf, when the latter is projected to the plane 
of these observables. This is a feature not uncommon in Bayesian analyses where the credibility 
regions map the 'volume' of scan points satisfying well enough a certain set of constraints, rather 
than representing iso-contours of the likelihood function, as is the case in frequentist analyses. The 
fact that the best-fit point is situated outside of the region of highest posterior probability for the 
considered observables simply tells us that, while for this point all the constraints are very well 
satisfied, it is also not very likely to obtain a similarly good value of the observable for similar 
choices of the input parameters. Particularly, in the case of Figs. 5|[a) and |^b)[ the majority 
of the points that satisfy all the constraints in the SC region present a Higgs mass in the range 
120-123 GeV, while for the other regions it is even less. 

Figure 2|[a)| also shows that the Higgs mass constraint favors the part of the AF region situated 
at m-i/2 <l 1500 GeV for a broad range of ttlq values. Besides, in Fig. ^ b) one can see that in the 
(^4.0; tan/3) plane the AF region spans a large range of positive Aq values and tan /3 is constrained to 
tan (3 > 40. In this region one loop corrections to the Higgs mass are driven up by a large stop mass 
("nrif ~ 2 — 4 TeV, as shown in Fig. a)) rather than large stop mixing which, given the allowed 



values of Aq, is minimal (see Fig. c [b)ty. Large values of MgusY push the Higgs mass close to the 
experimentally observed value, but not enough to reach it (for even larger MgusY) n-ai becomes too 
heavy and resonant annihilation of neutralinos is not efficient). This can be interpreted as a sign 
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Figure 5: Marginalized 2D posterior pdf in the (A, k) plane, for the CNMSSM constrained by the 
experiments listed in Table [T] for Case 1. The color code is the same as in Fig. [2j 



of some tension between the relic abundance and Higgs mass constraints in the AF region, which 
was investigated for the CMSSM in our previous paper |29j . The tension persists in the CNMSSM. 

Only a limited fraction of the FP/HB region (which appears as a narrow 95%-credibility 'tail' 
at mo 3> "1,1/2) survives, as was the case for the CMSSM, despite the fact that the relic density is 
well satisfied over there. As we pointed out in |29| . some tension with the 126-GeV Higgs boson 
mass arises not only in the AF but also in the FP/HB region where // e ff is smaller than anywhere 
else. In fact, small values of // e ff do not allow to obtain large negative At through running, as 
can be seen in Fig. ^[a)| where we show the posterior in the weak scale parameters At and fi e s 
(notice that At — Xt over all parameter space). For the chosen range of too, in the FP/HB region 
MgusY cannot be very large either, so that the correct Higgs mass cannot be reached. The region 
of high posterior thus moves up towards larger mi m , where the neutralino has still a non- negligible 
higgsino component, but -Mstjsy is large enough to give the correct Higgs mass. 

It is worth pointing out that Aq ~ is not realized in the CNMSSM. This is because the lightest 
pseudoscalar a\ becomes non-physical for such values. The mass of a\ is, for moderate and large 
values of tan/3, well approximated by ~ — 3ksA k |74| . A\ and A K are unified to Aq at the 
GUT scale, and A K barely runs, since the one-loop contribution to its /3-function is negligible. As 
a consequence, in the CNMSSM k and A K have always opposite signs and there are no points in 
the scan with Aq = or k = 0. 

We present in Fig. |^[b)| the marginalized posterior pdf in the (At, A T ) plane. The parameters 
show a clear linear correlation, which in the SC region (bottom left corner) results in large negative 
values for both observables, due to the fact that the correct Higgs boson mass requires large stop 
mixing, as discussed above. 

Figure [5] shows the 2D-posterior in the (A, k) plane. One can notice the known correlation 
between k and A [31], as \n\ cannot exceed a given value of A by too much without causing a 
Landau pole in RGE running and, on the other hand, A cannot exceed a given value of |k| by 
too much without lowering much below the observed value, where the likelihood becomes 
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Figure 6: Marginalized 2D posterior pdf in (a) the (tan/3, m^) plane and (b) the (m ai , tan/3) 
plane of the CNMSSM constrained by the experiments listed in Table [I] for Case 1. The outer 
contours enclose the 95% credibility regions and the inner contours the 68% credibility regions. A 
distribution of samples uniformly selected from our nested sampling chain is superimposed. 



negligible. In the SC region k is positive since A K is negative, while in the AF and FP/HB regions 
these parameters switch signs. Hence, the upper branch shown in Fig. [5] corresponds to the SC 
region, and the lower one to both the AF and the FP/HB regions. 

At 95% credibility, A and k show an upper bound, which changes over parameter space de- 
pending on the value of /U e g and loop corrections to m^, and again is due to the requirement 
of m/ ll ~ 126 GeV. When all the constraints are taken into account, the posterior shows that 
A < 0.16 for the upper branch, corresponding to the SC region, and A < 0.12 for the lower branch 
corresponding to the AF and the FP/HB regions. 

We want to re-emphasize that, when the global constraints are considered, Case 1 presents a 
very CMSSM-like character. The parameter A is small, and its effect on is insignificant, given 
that tan /3 > 4 over all of the regions of high posterior probability, as shown in Fig. 1 h) . In the 



SC region the lightest Higgs mass can assume values close to 126 GeV more easily than what was 
observed in our work on the CMSSM, but the reason lies in the updated value of the top mass used 
for the present analysis. 

In Fig. 6 ^a) we show the marginalized posterior in the (tan/3, m/ ll ) plane and in Fig. 6 [b) the 



posterior in the (m ai , tan/3) plane. We also overlap a distribution of samples uniformly selected 
from our nested sampling chain. Notice that the density of samples reflects their relative posterior 
probability. One can see in Fig. £ a) that, in the SC region (left island) the tension between the 



correct Higgs mass and the other constraints is much ameliorated. Figure 1 \b) shows confirmation 
of the CMSSM-like nature of this the pattern of high posterior mirrors the one found in 

Ref. 
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Figure 7: Marginalized 2D posterior pdf in (a) the (mo, m 1 / 2 ) plane and (b) the (Ao, tan/3) plane 
of the CNMSSM constrained by the experiments listed in Table [TJ for Case 2. The color code is the 
same as in Fig. [2j 



Case 2. The NMSSM allows more freedom in the Higgs sector than the MSSM, due to the 
extended number of parameters. Even in its partially constrained version, there is a possibility of 
obtaining a light h 2 , which we found to be a mixture of H u and S fields. A non-negligible singlet 
component creates the difference in the Higgs sector between this scenario and Case 1. In the rest 
of this section we analyze the consistency of an h 2 signal with the observed excess at the LHC. 



In Fig. "i a) we show the posterior pdf in the (mo, mi/2) plane for Case 2. The favored parameter 
space is now drastically reduced with respect to Case 1. Only the SC region, where the best-fit 
point is located, and the FP/HB region survive the requirement of having mh 2 ~ 126 GeV. 

In Fig. I'fflb)! we show the posterior pdf in the (Ao, tan/3) plane. In both regions the range 
allowed for Aq does not extend much from zero, -400 GeV < Ao < 400 GeV0 The reason is that 
\A K \ is now limited to values less than 400 GeV by our requirement on the mass of h 2 . Specifically, 
the masses of the two lightest CP-even Higgs bosons at the tree level are approximately given 

byHU 

m 2 hl 2 « i |m| + A(ks) 2 + ksA k T s] [Mf - 4( KS ) 2 - ksA k ] 2 + A-term j , (17) 
where the A-term under square root includes the usual correction to the tree level masses and 



depends on the other parameters of the Higgs sector. It is clear from Eq. (17), that in the regime 
where \k\s < Mz, the mass of h 2 is of order Mz, while m^ scales as \k\s. When approaches 
zero, the A-term is C(< 10 6 GeV), while the first term under square root is roughly of the order of 
M| ~ 10 8 GeV or larger. The A-term can therefore be neglected, and the requirement of having 
m \ 1 — translates into the approximate relation \A K \ < 4|k|s. On the other hand, in the regions 
where \k\s > Mz, rrih 2 ~ us and ~ Mz- Values of m^ 2 much greater than 126 GeV are 



Note that a gap around Ao ~ in Fig. "1 'b) comes from the physicality condition for an,, as discussed above. 
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Figure 8: (a) Marginalized 2D posterior pdf in the (m^, mh 2 ) plane of the CNMSSM constrained 
by the experiments listed in Table [T] for Case 1 . |(b)| Marginalized posterior pdf in the (m^ , rrih 2 ) 
plane of the CNMSSM for Case 2. The color code is the same as in Fig. [2j 



forbidden by the likelihood function, so that \k\s presents an upper bound, which translates into 
an upper bound on \A K \. 

Since in most of the parameter space s is very large, and A and k are correlated, the scan also 
shows upper bounds for A and k. In the SC region A is very small, A < 0.01, while in the FP/HB 
region it can assume slightly larger values, A < 0.04. Notice that the upper bound on \k\s does not 
depend on the other input parameters, but the bound on n (or A) does, particularly on the value 
of //eg. In the SC region ix e s > 600 GeV while in the FP/HB region /U e g ~ 200 GeV. 

Since in the SC region A and k are bounded to be much smaller than in Case 1, the lightest 
neutralino there is a nearly pure singlino and very light (the neutralino mass matrix with the 
convention used in this paper can be found, e.g., in [31]). The corresponding tan/3 assumes larger 
values (tan j3 ~ 30 — 35, favored also by other constraints, e.g. 5 (g — 2) ) than in the same region of 
Case 1, or of the CMSSM. For smaller values of tan/3, on the other hand, the lighter stau becomes 
the LSP. 

In Figs. £|[a) and £[b) we show the 2D posterior in the (m^, mh 2 ) plane for Case 1 and Case 2, 
respectively. One can see that in Case 1 m^ 2 3> m^, while in Case 2 the 68% credibility region 
shows a preference for 90 GeV < < 115 GeV. 



Since in Case 2 mh 2 ~ 126 GeV, 



in 



O-l 



assumes values quite small in the favored regions, 



■m, 



«i 



< 



500 GeV. While it might appear that such a case is excluded by the recent H — )■ tt searches at 
the LHC [75], we point out that the limit does not apply to this case, as the relevant couplings for 
a\ — > tt and a\ — > bb are suppressed by the pseudoscalar mixing angle Op, i.e., when | cos#p| < 
1 [76]. We have checked that | cos#p| is very close to zero over the regions of interest. 
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Figure 9: |(a)| Marginalized 2D posterior pdf in the (R hl (jj), R hl (ZZ)) plane of the CNMSSM 
constrained by the experiments listed in Table [I] for Case 1. (b) Marginalized posterior pdf in the 
(i?/i 2 (77), Rh 2 {ZZ)) plane of the CNMSSM for Case 2. The color code is the same as in Fig. [2j 



Case 3 (two degenerate light Higgs bosons) is in fact a subset of Case 2 when it comes to the 
allowed parameter space, which is uniquely determined by the requirement of rrih 2 ~ 126 GeV as 
was explained above. The posterior pdf in this case is very similar to the one showed in Figs. 



and and the best-fit point is again situated in the SC region. However, this case presents also 
some characteristic features: 

• h\ is a mixture of and H u , with a predominance of the latter, while /12 is mainly singlet-like, 
with a small fraction of H^. 

• Both A and k are now much more limited and very close to zero. The reason for this is simple: 
in order to have the two lightest CP-even Higgs bosons almost degenerated in mass, one needs to 



minimize the difference m\ 2 — from Eq. (17). This yields negligible values for the parameter 
k and, consequently, also for A. 

4.3 Impact of the cross section rates 



In Fig. S|[a) we show the 2D posterior pdf in the (R^iZZ), ^^77)) plane for Case 1. As one can 



see, the Higgs boson in this case is SM-like. Actually, it is known [30] that the enhancement of the 
signal strength in the 77 decay channel observed by both CMS and ATLAS cannot be obtained for 
the values of A that are favored by the global scan. As a consequence, 12^(77) cannot be fitted 
perfectly in the CNMSSM. This discrepancy effectively adds two units of x 2 homogeneously over 
the allowed parameter space, but it does not alter the posterior distribution. For the same reason 
Rh^iZZ) cannot be perfectly fitted either, though its contribution to the total x 2 is smaller than 
0.5 units of x 2 , making this observable equally ineffective to constraining the posterior. 

In Fig. £|[b) we present the posterior distribution for Case 2. Once again, Rh 2 (ll) can hardly 



become larger than one over the allowed parameter space. The 95% credible region lies far from 
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the central value of the observed enhancement and, in fact, even covers values lower than in Case 1. 
Rh 2 (ZZ) presents similar behaviour, although the suppression of the reduced cross section is highly 
welcome for this observable, as it places the calculated value closer to the rate observed at CMS. 
Smaller than one signal rates indicate less of a SM-like character for h,2, which is caused by the the 
suppression of the SM-couplings induced by its increased singlet component. 

The posterior distributions presented in Figs. l [a)| and ^b)| indicate that, in both Case 1 and 



Case 2 it is in general extremely difficult to obtain the signal enhancement in the 77 channel. The 
scan naturally tends to stay in the regions of parameter space favored by all constraints. It is 
therefore no surprise that among the points scanned for Case 1 only two presented a 77 rate in the 
range 1.2 — 2, thanks to the reduced coupling of the signal Higgs boson to the bottom quarks. Such 
points present x 2 contributions to the relic density of order several tens, and the x 2 contribution to 
BR (B s — > fi + /j>~) is of order 100. In Case 2 we found a dozen such points, for which the contribution 
to the relic density is even worse. 

In Case 3 one could expect to obtain an enhancement of Rh si (77) by adding the individual 
rates for both almost degenerate light scalars. However, the posterior pdf in the (Rh 1 +h 2 (,1l) j 
Rhi+h 2 (ZZ)) plane is remarkably similar to the one shown in Fig. £ [a)| for Case 1, due to the large 



singlet component of /i2- In fact, in Case 3 we were not able to find a single point with the enhanced 
77 rate. Since Case 3 is a subset of Case 2 in terms of the favored parameter space, and the rates 
in the 77 and ZZ channel do not show interesting features, we will not consider it any further. 

4.4 Prospects for DM direct detection and BR(B S — > ) 

In this subsection we will discuss the impact of limits from direct DM searches on the allowed 
parameter space of the CNMSSM. This kind of experiments is complementary to direct LHC 
SUSY searches, as they are capable of testing neutralino WIMP mass ranges beyond the current 
and future reach of the LHC, and therefore could add new pieces of information to the global 
picture. 

At present the most stringent limit on the spin-independent cross section cr^ 1 comes from 
XENON100 [77J. In supersymmetric models it can then be plotted as a function of the neu- 
tralino mass in the form of an exclusion limit in the (m x , dp 1 ) plane. In this study, we did not 
include the upper bound on ct^ from XENON100 in the likelihood function, since it was shown 
in |41j that, given the large theoretical uncertainties (of order 5 to 10), it typically does not affect 
the shape of the posterior pdf. 



In Fig. 1C ^a) we present the posterior pdf in the (m x , (Tp 1 ) plane for the Case 1. The solid 
red line shows the most recent 90% CL upper exclusion limit by the XENON100 Collaboration. 
One can see that the two high posterior probability regions (from SC on the left and AF on the 
right) are located well below it. A small 95% credible region (corresponding to the FP/HB region 
with the lightest neutralino being a mixture of bino and higgsino) is almost (but not yet entirely) 
excluded, even taking into account significant theoretical uncertainties. One can also notice that 
XENON1T, a future ton-size DM detector (projected sensitivity represented as the dashed gray 
line), will be in a position to test the 68% credible regions. 



The impact of the XENON100 limit in Case 2 is even stronger. In Fig. [Kp)l we show the 



posterior in the {m x , a^ 1 ) plane. The high probability region above the experimental line presents 
the features of the FP/HB region: /x e g ~ 200 GeV, and the lightest neutralino is a mixture of bino 
and higgsino. As a consequence, the spin-independent neutralino-proton cross section tends to be 
larger, so that the region has the potential of being ruled out soon. Since we force m^ 2 ~ 126 GeV, 
the favored region of m x is pushed down to < 240 GeV. This is important as it implies that a 
significant fraction of the parameter space for the case in which the role of the SM-like Higgs is 
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Figure 10: Marginalized 2D posterior pdf in the (m x , a^ 1 ) plane of the CNMSSM constrained 



by the experiments listed in Table [T] in (a) Case 1 and (b) Case 2. The solid red line shows the 
90% CL exclusion bound by XENON100 (not included in the likelihood), and the dashed gray line 
the projected sensitivity for XENON1T. The color code is the same as in Fig. [2j 



played by the second-to-the-lightest CP-even scalar is strongly disfavored. 

Notice that the SC region is not shown in Fig. 1C [b) since the neutralino is a nearly pure 
singlino there, so that c^ 1 is several orders of magnitude below the XENON100 bound (and below 
the range shown in the figure). However, we would like to point out that in Case 2 (and Case 3) 
strong constraints on the full parameter space can be placed as a result of the interplay between 
the limits provided by two completely different experiments - that test different observables by 
means of different experimental techniques - namely LHCb and XENON100. This is illustrated in 
b) where we show the posterior pdf in the (BR(B S — > /x + /i _ ) , Tp 1 ) plane for 
respectively. The solid blue horizontal line shows the minimum 90% CL upper 
bound on a^ 1 by XENON100 (obtained at m x ~ 50 GeV), and the dotted blue horizontal line the 
corresponding projected sensitivity for XENON1T. The pink vertical band encompasses the lcr 
experimental uncertainty on the recent LHCb measurement of BR(B S —> [35]. Figure 11 'b) 

shows that, for Cases 2 and 3, in the SC region BR(B S — > /i + /i~) is strongly enhanced, due to the 
large values assumed by tan/3 there, and it could be excluded by the next updated results from 
LHCb. 



Figs. 



11 



11 



and 

Case 1 and Case 2 



4.5 Fine-tuning 

In this subsection we will address the issue of fine-tuning. Note that we will not delve into it, 
nor will we discuss which values of fine-tuning are acceptable or not from the point of view of 
naturalness. Our aim here is to simply present an estimate of fine-tuning (provided as an output 
by NMSSMTools v3.2.1) for the allowed parameter space of the model, leaving aside the discussion 
of the viability of the model itself, which would be a matter of personal prejudices. 

The mass of the Z boson (which determines the EW symmetry breaking scale) can be expressed 
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Figure 11: Marginalized 2D posterior pdf in the (BR (B s -)• ^ + n~), a^ 1 ) plane of the CNMSSM 
constrained by the experiments listed in Table [T] in | (a) | Case 1 and |(b)| Case 2. The color code is the 
same as in Fig. [2 The solid blue horizontal line shows the minimum 90% CL upper bound on cr^ 1 by 
XENON100 (not included in the likelihood), and the dotted blue horizontal line the corresponding 
projected sensitivity for XENON1T. The pink vertical band shows the la experimental uncertainty 
on the recent measurement of BR (B s — > [35] (which is included in the likelihood). 



in terms of the supersymmetric parameters through the minimization condition of the Higgs po- 
tential, 



Ml _ _ 2 rn 2 Hd - m 2 Hu tan 2 p 
2 M + ~ (tan 2 /3-l) 



P + -^7 20 ^ ■ ( 18 ) 



The fine-tuning problem of the MSSM [78] amounts to the fact that the parameters m,H u , mH d 
and fi need to be simultaneously tuned to a high precision to reproduce the correct value of Mz- 
In the NMSSM fi is replaced by fi eS . 

In the framework of the unified theory, at the GUT scale mn u = m H d = mo- They are then 



evolved down to the EW scale by means of the RGE's. Therefore, the right-hand side of Eq. ( 18 ) 
depends on the parameters of the model at the GUT scale. 

The measure of fine-tuning associated with the parameter pi of the model is defined as [79l [7S] 



Pi 



dlogMf 



dlogp 2 



(19) 



In the CNMSSM, pi = {mo, mi/ 2 ,Ao, A, k}. We do not present the fine-tuning measure due to the 
top-quark Yukawa coupling, following the approach adopted in [73] , and quantifying the fine-tuning 
associated only with the supersymmetric parameters. We define the fine-tuning measure A for a 
given model point as the maximal contribution to fine-tuning among all the model's parameters 
for that point, 

A = Max(A Pi ). (20) 
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Figure 12: Iso-contours of the fine-tuning measure A (dashed red lines) superimposed to the 2D 
posterior pdf for (a) Case 1 and (b) Case 2. The color code is the same as in Fig. [2j 



In Figs. 12;|^a) and 12j[b) we present the iso-contours of the fine-tuning measure A in the (m-o, 
m i/2) plane for Cases 1 and 2, respectively. The iso-contours reflect the value of A for the vast 
majority of the points included. They are superimposed on the 2D posterior distributions. Case 1 
presents very CMSSM-like behaviour, where smaller fine-tuning can be achieved only in the FP/HB 
region due to the relatively low values of // e ff. Note that for the same reason A is larger in the SC 
region which is characterized by larger /x e ff. On the other hand, Case 2 is less fine-tuned (A < 50 
in a vast region of the parameter space allowed at 2cr), which is a reflection of the fact that the 
parameter space is already highly constrained by the requirement of mh 2 ~ 126 GeV. In Figs. l^[a) 
and l^[b) we show for what percentage of the total number of allowed points each of the p^s yields 
maximal fine-tuning. For example, k gives maximal fine-tuning for 94% of the points in Case 1, 
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Figure 14: Marginalized 2D posterior pdf in (a) the (mo, m 1( / 2 ) plane and (b) the (Ao, tan/3) plane 



of the CNMSSM constrained by the experiments listed in Table [TJ for Case 1 without (g — 2)^ and 
with /j, e g < 0. The color code is the same as in Fig. [2j 



while Aq is the main contributor for the vast majority of the points in Case 2. 
4.6 Relaxing (g — 2) M and the case of negative /j, e g 

In our recent study of the CMSSM [29] we considered the effect of relaxing the 5 (g — 2)^ con- 
straint. The reason was based on the observation that the poor fit of the CMSSM was the result 
of basically that single constraint, which simply could not be reproduced after including especially 
direct superpartner mass limits from the LHC. This appears to be a general feature of all SUSY 
models where slepton and squark masses are assumed to be comparable, including simple unified 
models like the CMSSM or CNMSSM. We simply state it as a conclusion reached by many studies. 
While we do not feel to be in a position to comment on the reliability of theoretical calculations 
of the SM (g — 2)^ which are strongly affected by non-perturbative effects related to low-energy 
strong interactions, especially the hadronic light-by-light contribution (for more details see, e.g., 
the Introduction of |29j and references therein), we believe that it thus makes sense to consider 
global scans where the constraint is removed. We showed in [29J that for the CMSSM the better 
fit was obtained with fi < thanks to the constraints from r>physics, which present a much better 
X 2 for this choice of parameters. 

Relaxing the (g — 2)^ constraint in the CNMSSM while keeping // e fF positive has no apparent 
impact on the posterior distributions both for the parameters of the model as for the measured 
observables. Such a behavior was to be expected and has already been observed in [29]. We also 
confirm that the better overall fit is obtained with // e g < for basically the same reasons as in the 
CMSSM. 

we show the 2D posterior pdf in the (mo, mi/v) plane for Case 1, where we 
ignored the constraint from 5 {g — 2) and set sgn(^ e g) = —1. Differently from the CMSSM, where 



In Fig. 14' a) 
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Figure 15: Marginalized 2D posterior pdf in (a) the (m 

(Tp) plane of the CNMSSM constrained by the experiments listed in Table 

(9 



plane and (b) the (BR (B s — > ), 
for Case 1 without 



2)„ and with /j, e g < 0. The color code is the same as in Figs. 10 and 



able |il 
d|ll 



one could observed a clear predominance of the AF region, in the CNMSSM the SC, AF and 
FP/HB regions of parameter space now seem to be equally probable. This is not likely to be an 
intrinsic difference between the CMSSM and the CNMSSM; it is more probably due to the fact 
that the constraints are implemented with different numerical tools in the two models, as explained 
exhaustively at the beginning of Sec. |4| We note here that the global likelihood becomes flat over 
the regions of parameter space allowed by the relic density, and very little information can be 
extracted from the posterior. 

In Fig. l4|[b) we show the same posterior pdf in the (Aq, tan/3) plane for Case 1. As was the 
case for the CMSSM, one can see that the distribution of tan /3 tends to favor slightly smaller values 
than in the positive /U e ff case, particularly in the AF region. 

11 a) we show the 2D posterior pdf in the (m x , ap^ 1 ) plane for Case 1. Since the posterior 
HB region is much more extended with respect to the positive case, a large region 
of parameter space lies above the XENON100 bound, and has the potential of being tested with 
modest improvements in sensitivity. On the other hand, as was the case for the CMSSM, the AF 



In Fig. 
in the FP, 



region (m x > 



400 GeV and erf 



< 10 pb) is not likely to be further constrained by the new 
spin-independent cross section measurements planned for the next years, including XENON1T. 

We shall analyze the x 2 contributions from the individual constraints in the next subsection. 
Here we only repeat that, when the (g — 2)^ constraint is relaxed, the better \ 2 can be obtained 
for /i c g < thanks to the improved fit to BR(B S — > and BR (B — > X s 7). To illustrate this 

feature we show in Fig. II: b) the 2D posterior pdf in the (BR(B S — > o^ 1 ) plane for Case 1. 

When /i e fj < 0, BR (B s — > [i + receives negative SUSY contributions that improve the fit over all 
parameter space, and particularly in the AF region, where the calculated values are pushed below 
the SM value. 
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4.7 Mass spectrum and the best-fit points 



In Figs. It ^a) and It \b) we show the ID marginalized posterior pdf's for the SUSY mass spectrum 
in Case 1 and Case 2, respectively. The narrow lines indicate the 95% credibility regions, and the 
thick bars the 68% credibility regions. 

Figure It [ a) | shows the CMSSM-like character of Case 1, since /13, 02 and %5 are a h heavy and 
effectively decoupled from the low scale spectrum. For Case 2, Fig. It \b) shows that requiring 
rrih 2 ~ 126 GeV forces the particles of the Higgs sector to be quite light, while the sfermions remain 
heavy, although they tend to be lighter than in Case 1, since is lighter and it requires smaller 
loop corrections. The latter affect the scale of the stop and, through the assumption of unification, 
the other sfermions. 



In Fig. It c) we show the mass spectrum for Case 1 where we neglect 5 (g — 2)^ and set ^ e ff < 0. 
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Figure 16: SUSY mass spectra of the CNMSSM with ^ eff > for 



Case 1 and 



(b) 



Case 2. 



(c) SUSY spectrum with < and no {g — 2)„ in Case 1. The narrow lines indicate the 



credibility ranges and the thick bars the 68% credibility regions. 
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Figure 17: Individual x 2 contributions to the best-fit points of the scans considered. 



Again we confirm the CMSSM-like nature. 

In Fig. 17 we show the individual x 2 contributions of the applied constraints to the best-fit 



point, for the cases considered in this study. As could be expected, in Case 1 the main contribution 
comes from 5 {g — 2) , once more confirming the CMSSM-like character of this scenario. In Case 
2 the contribution from 5 (g — 2) is comparable to the one from BR (B s — > due to the 



singlino character of neutralino which favors higher tan/3, as explained in Sec. 4.4 When [i e R > 0, 



relaxing the (g — 2) M constraint does not change significantly the contributions to Xmin due to the 
other constraints. This shows that the x 2 due to (g — 2) M is fairly homogeneous throughout the 
parameter space not yet excluded by direct LHC SUSY searches. However, for negative /i e fi the 
total x 2 improves by around 1-2 units, thanks to the improved fit to BR (B — > X s 7), for which 
the contribution from the chargino-stop loop changes sign and tends to enhance the branching 
ratio, towards the experimental value |80| . As a consequence, the overall fit to the experimental 
measurement improves. 

One can finally notice that the contributions to Xmm due to Rh Bi (77) are sizeable in all three 



cases, as was explained in Sec. 4.3 



In Table [3] we show the parameters of the best-fit points obtained in our scans. 



5 Summary and conclusions 

In this paper we have presented the first global analysis of the CNMSSM which included the 
measurement of the mass and decay cross sections of the Higgs-like resonance observed at the LHC, 
the first evidence of nonzero BR (B s — > at LHCb, SUSY mass limits from direct searches, 

as well as the updated value of the top mass. We discussed in detail the possibility of either one 
of the lightest CP-even Higgs bosons of the model playing the role of the SM-like Higgs (Case 1 
and Case 2 in the text) and that of the observed excess being due to a combination of two of them, 
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tan /3 
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m h Bis 


A 




With (g - 2) M 


Case 1, ^ e ff > 


340.9 


911.9 


-2546 


12.0 


0.017 


124.5 


427 


17.91 


Case 2, /j, c s > 


102.8 


803.6 


-295.0 


32.1 


0.006 


123.5 


396 


18.45 


Case 3, /j, c s > 


124.7 


958.8 


-356.7 


33.7 


0.001 


120.3, 126.2 


480 


20.08 


No (g - 2)„ 


Case 1, /U e ff < 


410.2 


913.3 


-2764 


14.8 


0.017 


123.9 


455 


6.40 



Table 3: CNMSSM input parameters, Higgs masses and fine-tuning measure for the best-fit points 
of our scans. Masses and Aq are in GeV. 

both lying within the range ~ 123 — 129 GeV, after inclusion of a 3-GeV theoretical error on the 
Higgs mass calculation (Case 3 in the text). 

Contrary to what is generally believed, we found that the CNMSSM does allow a SM-like 
h\ as heavy as 125 GeV, especially in the stau coannihilation region where the scan obtains a 
good compromise between the relic density and the h\ mass constraint owing to large negative Aq 
values. The calculated value of m/ ll is larger than what was observed for the CMSSM in the same 
region thanks to the increased value of the top pole mass. The overall fit is somewhat spoiled by 
the requirement of having (77) consistent with the value reported by the CMS collaboration, 
since it is virtually impossible to obtain an enhancement in the cross section rate in Case 1. A 
similar conclusion can be drawn for Case 2 where, in turn, h 2 is required to have a mass around 
126 GeV, which can also be achieved but at the cost of a poor fit to other observables. Case 3, 
with almost degenerate h\ and h 2 was found not to be as interesting as anticipated, since the 
combined ^^+^(77) almost never exceeds 1, besides the fact that the high credibility posterior 
regions mimic those of Case 2. 

We also provided estimates of fine-tuning due to the various input parameters of the model in 
the form of iso-contours in the (mo, ^1/2) plane for Cases 1 and 2. We noted that the maximum 
fine-tuning for most of the parameter space comes from two different sources in the two cases. 

Finally, we assessed the effects of abandoning the (g — 2) M constraint since it cannot be repro- 
duced in the CNMSSM, and more generally SUSY models with slepton-squark unification. In this 
case the overall fit actually improves considerably for fi c g < 0, due to a better agreement of the 
model's predictions for 6-physics observables with experimental data, similarly to the CMSSM. 
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